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Abstract; [ Aim] Apis cerana cerana is an important germplasm resource and pollinator. This study aims to 


explore the genetic diversity and population structure of A. cerana cerana in the Yimeng Mountains, Shandong 


province, East China so as to provide a theoretical foundation for the utilization and conservation of germplasm 


resource for A. cerana cerana. [Methods] In this study, 114 broods of A. cerana cerana from 7 localities in the 


Yimeng Mountains were selected and 36 morphological indexes were measured according to Ruttner' s analytical 


method. And the genetic diversity of A. cerana cerana in the Yimeng Mountains was evaluated with 11 


microsatellite loci. [Results] The results of morphological analysis showed that the average body length and the 


average forewing length of A. cerana cerana foraging workers in the Yimeng Mountains were 12. 064 — 13.351 


mm and 8. 198 — 8. 694 mm, respectively. Fifty-eight alleles were found in the 11 microsatellite loci. The 


number of alleles per locus ranged from 2 to 17. The average expected heterozygosity and polymorphism 
information content (РІС) of all loci were 0. 3115 and 0.2872, respectively. The average number of alleles in 
all the populations ranged from 2. 4545 (ST-AQ) to 4.0000 ( BHY) , and the average expected heterozygosity 
ranged from 0. 1916 to 0. 3397. The average genetic differentiation measured as an Fst value was 0. 048. The 
Nei's genetic distance ranged from 0. 0092 ( XL-DLZ to XZ-XLZ) to 0. 1000 ( XL-DLZ to XL-DJW ). 


Besides, the cluster analysis results showed that all the populations in the Yimeng Mountains could be grouped 


into three categories. [Conclusion] The population structure of A. cerana cerana populations in the Yimeng 


Mountains is similar. The results suggest that the gene flow occurs in different populations. Compared with the 


populations in other positions, A. 
morphology and genetic structure. 


cerana cerana in the Yimeng Mountains shows certain specificity in 


Key words: Apis cerana cerana; morphology; microsatellite; genetic diversity; Yimeng Mountains 


1 INTRODUCTION 


The Chinese honeybee, Apis cerana cerana, is 
the geographic subspecies of Apis cerana. Compared 
with the Western honeybee ( Apis mellifera) , it has 
an acute sense of smell and a strong resistance to 
mites and can forage a wide range of nectars and 
pollens ( Li et al., 2008). A. cerana cerana has a 
wide distribution in China. However, the qualities of 
A. cerana cerana in China have decreased to less 
than 10 000 broods 


competition since the introduction of A. mellifera. 


because of interspecies 


The Yimeng Mountains in Linyi city, Shandong 
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province, East China, has a temperate continental 
monsoon climate. This area has sufficient sunshine 
and the weather is relatively moderate with an annual 
average temperature of 13. 1. Тһе Yimeng 
Mountains has a long history of А. cerana cerana 
rearing with traditional farming methods, such as vat 
farming. 

There are many research methods for examining 
honeybee genetic diversity. Different methods can 
reveal genetic information from different angles and 
levels. A morphological method can describe the 
genetic variation by measuring morphological and 
phenotypic characteristics ( Ruttner, 1988; Li and 
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Wang, 2004). The morphological characteristics 
mainly include body length, length of femur, cubital 
index, etc. All the morphological characteristics are 
important for line identification. Before the 1960s, 
morphological analysis was an important method for 
variety identification and genetic structure analysis in 
honeybee. Morphological analysis is a common tool 
for determining honeybee genetic diversity ( Ruttner, 
1988). 

Due to their high polymorphism, 
specificity, abundance and random distribution over 
the genome апа inheritance, 
microsatellites are most commonly used to assess 


locus 
codominant 


population structure and diversity ( Weigend and 
2001 ). The Honeybee 
Sequencing Consortium reported that there are more 
than 2 000 microsatellite markers. The average 
distance between microsatellite markers is 2. 1 cM 


Romanov, Genome 


( The Honeybee Genome Sequencing Consortium, 
2006 ). According to FAO 


determining classic genetic distances using neutral, 


recommendations , 


highly polymorphic microsatellite markers is the 
method of 
relationships and breed differentiation ( Crooijmans 
et al., 1996). This methodology also provides 


choice for investigating genetic 


information for establishing preservation priorities for 
livestock breeds ( Barker, 1999). Estoup et al. 
(1993) analysed plesiomorphic bee species based on 
two microsatellite markers. Franck et al. ( 2000) 
proved that the honeybee in the Near East may be 
the fourth bloodline of A. mellifera. 

Analysing the genetic diversity of A. cerana 
cerana in the Yimeng Mountains could provide a 
theoretical foundation for the research on the 
conservation and reasonable utilization of honeybee 
populations. In this study, essential data of 11 
microsatellite markers and 36 morphological indexes 
from seven locations in Yimeng Mountains are 
provided. The results may be useful to understand 
the genetic differentiation in this important genetic 
efficient 


resource and contribute to more 


conservation. 


2 MATERIALS AND METHODS 


2.1 Biological specimens 

We chose the localities based on the honeybee 
breeding method, which was vat breeding or soil 
breeding, and all the bees were collected from the 
field. We chose seven localities, and their detailed 
information was shown in Table 1. 


Table 1 Sampling information of Apis cerana cerana in the Yimeng Mountains, Shandong 





Populations Sampling locality за Қ” pow oru 
XZ-XLZ Xiaoliuzhuang Village, Xuezhuang Town, Linyi City 2017/7/7 35?21'36"N, 11821026 Е 200 20 
DTZ-ZJZ  Zhoujiazhuang Village, Datianzhuang Town, Linyi City 2017/1/1 35?27'38"N, 118?0'45"E 350 15 

BHY Baihuayu, Linyi City 2017/7/7 35°33'31"N, 117555'35"E 550 20 
XL-DJW Dajiawan Village, Xili Town, Zibo City 2017/71/29 36°0'35"N, 118?14'58"E 338 15 
XL-DLZ Daliuzhuang Village, ХШ Town, Zibo City 2017/71/29 36°0'12"N, 118?15'14"E 385 15 
ZZ-ZHZ Zhonghanzhuang Village, Zhongzhuang Town, Zibo City 2017/1/29 36?2'26"N, 118?17'9"E 316 10 
ST-AQ Anging Village, Sitou Town, Weifang City 2017/71/29 36?18'43"N, 118?31'32"E 174 19 


To obtain similar aged bees, we collected 
samples at the door of the beehive. In total, 114 
broods from seven localities in the Yimeng Mountains 
were collected. Ten foraging workers were collected 
from each brood, and the samples were placed in 
75% ethyl alcohol as materials for the following 
morphological study. Moreover, one more worker 
bee was collected from each brood and placed in 
95% ethyl alcohol as materials for the following 
microsatellite study. DNA was extracted according to 
reports by Ji et al. (2009) and then preserved at 
– 20% after the content and purity were tested. 
2.2 Morphological index measurement and 
analysis 

Ten bees were selected randomly from each 
brood. Then, all the bees were dissected and their 


36 morphological indexes, including body length, 
length of femur, cubital index, etc., were measured. 


The measurement standards referred to Wang 
(2007). 
2.-3 Microsatellite primer selection and 


multiplex PCR amplification 

According to the results of Ji et al. (2014) , 11 
pairs of microsatellite primers were selected ( Table 
2). The 11 pairs of microsatellite primers were 
labelled with FAM (blue), HEX (green) or 
TAMRA ( yellow) fluorochromes at the 5'-terminus. 
Then, the primers were combined according to their 
base length ( Ji et al., 2014). The primer 
combination information was included in Table 3. AII 
the primers were synthesized by Sangon Biotech 
( Shanghai). PCR products were obtained in а 20 uL 
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Table 2 Microsatellite primers used in this study 
Locus Primers Primer sequences (57-37) Repetitive unit GenBank accession no. 
K 29-CL1229. Contig33F CACAAGGAGTTCTTTCTCGACAT (AT)8 NW_003377929. 1 
29-CL1229. Contig33R TTCCACCGTATTGGAATTGATTA 
L 33-CL1360. Contigl4F AGTTCGACAGACCAAGCTGTAAG (TG)6 NW_003378041. 1 
33-CL1360. ContigI4R CTTGTGCGAACATAGAACAAAGA 
M 1-CL1229. Contig27F ATTCTACGATACGGCACGTTACA (AT)8 NW_003378122. 1 
1-CL1229. Contig27R ATTGTGAGTGGGTACAGTTTCGT 
P 7-CL1278. Contig] F TGTGATAGAATCGTGTGTCTTACG (TG)7 NW_003378123. 1 
7-CL1278. Contig] R GCCAATGTAATAACGATAGTGCC 
Q 11-CL133. Contig41F AGAGGGGGAAGAAAAGAAAAGAA (TCG )6 NW_003377943. 1 
11-CL133. Contig41 R GACCGGTTAATCACTGGTGTC 
R 14-CL1549. Contig3 F CGTCCAAGTCTGTACTCCAGC (CCG)5 NW. 003378143.1 
14-CL1549. Contig3R AGTAAGTGACGAACACCGTATCG 
5 17-CL1462. Contig5 F AACAAATAGGCCTCGTTAATGTTT (AC)9 NW. 003377909. 1 
17-CL1462. Contig5R CGGCTGCTTCTTTCTCTTTTTAT 
T 18-CL1470. Contig3 F AATCGATAGACGATGAATTTGGA (ТАА)5 NW. 003377909. 1 
18-CL1470. Contig3R GGACGAAATGGATCAAAGTAAGA 
X 39-CL308. Contig6F TAGAGGTGGAAGGCAATTATCTG ( AT)6 NW. 003377999.1 
39-CL308. Contig6R AAATTATGTGATTATTACAATGAGTACC 
W 37-CL1293. Contig] F GATAGTTACACGCCGCTTAAATG (AT)8 NW_003377999. 1 
37-CL1293. Contig] К GATTACTAGATCCAGCGCACATC 
y 51-CL1114. Contigl5F TCGAATGTTAGGATTTTCGTGAT (CG)6 NW_003378039. 1 
51-CL1114. Contigl 5R ATGCGAGTACGTGTGAATGTATG 
Table 3 Microsatellite primer combinations for multiplex PCR in this study 
Dus Ris быры с ПГУ” 
accession no. (bp) (©) position used 
P 7-CL1278. Contigl LOC410851 NW_003378123. 1 FAM (blue) 159 59.5 LG2 
y 51-CL1114. Contigl Mblk-1 NW_003378039. 1 HEX ( green) 145 59.9 1615 P&y&X 
X 39-CL308. Contig6 LOC412784 NW 003377999.1 TAMRA (yellow) 114 57:9 LG9 
5 17-CL1462. Contig5 complexin NW_003377909. 1 FAM (blue) 159 59.9 LG4 
T 18-CL1470. Contig3 10С413936 NW_003377909. 1 HEX ( green) 142 59.6 LG4 S&T&W 
W 37-CL1293. Contigl LOC100749790 | NW 003377999.1 TAMRA (yellow) 89 60.1 LG9 
K 29-CL1229. Contig33 — LOC100870946 | NW 003377929.1 FAM (blue) 156 60.1 LG7 K 
L 33-CL1360. Contigl4 PMCA NW. 003378041.1 FAM (blue) 151 59.7 LG8 L&M 
M 1-CL1229. Contig27 LOC100870946 =NW_003378122. 1 HEX ( green) 137 60.3 ІСІ 
R 14-CL1549. Contig3 КСА NW. 003378143.1 FAM (blue) 150 60.2 LG3 R&Q 
Q 11-CL133. Contig41 LOC100577278 NW _003377943. 1 HEX ( green) 136 60.1 LG3 





reaction volume containing 10 x Buffer 2. 0 pL, 
MgCl, (25 mmol/L) 1.0 uL, dNTPs (10 mmol/L) 
0.5 uL, forward primer (10 pmol/uL) 1.0 uL, 
reverse primer (10 pmol/wL) 1.0 uL, Taq DNA 
polymerase (5 U/uL) 0.2 uL, worker bee DNA 
template (1.0 ng/wL) 1.0 uL and ultrapure water 
13.3 рі. PCR was performed according to the 
following program: (1) 5 min at 95%; (2) 30 
cycles of 50 s at 95*C , 50 в at 59 – 62%, and 50 s 
at 72°C ; and (3) 10 min at 72%. Six microliters of 


the products from the last step was used to perform 
agarose gel electrophoresis for 20 min at 120 V to 
examine the PCR products. Short tandem repeat 
( STR ) genotyping of the PCR products was 
performed by Sangon Biotech and the atlas files were 
generated by GeneMapper 4. 0. 
2.4 Statistical analysis 

Discriminant analysis and clustering analysis 


were performed based on the morphological indexes 
with SPSS 22.0. The allele frequency, and observed 
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heterozygosity (Ho) and expected heterozygosity 
(He) (Nei, 1987) for each population across the 
loci and for each locus across the populations were 
estimated with the Microsatellite-Toolkit for Excel 
based on the results of microsatellite genotyping. 
The polymorphism information content ( PIC) for 
each locus was obtained according to Botstein et al. 
(1980). The F-statistics indexes ( Wright, 1978) 
were estimated in the form of Fit, Fst and Fis, the 
sample-based, respective these 
parameters proposed by Weir and Cockerham 
(1984), as implemented in the FSTAT program. 
The significance of the F-statistics was determined 
from permutation tests with the sequential Bonferroni 
procedure applied over the loci ( Goudet, 2002). 
The Fst values among the populations and effective 
individual migration in each generation were 
calculated with the GENEPOP program ( Raymond 
and Rousset, 1995). The Nei' s genetic distance 
and genetic identity were calculated with Popgene 
32. The phylogenetic trees were created based on 
Nei’ s genetic distances using the neighbour-joining 
method. The genetic structure of the population was 
estimated with the Structure program and the optimal 
cluster numbers were calculated using Structure 
Harvester (http: // taylor0. biology. ucla. edu/ 
structureHarvester/). The genetic structure figure 


estimators of 


was drawn with Distructure. 


3 RESULTS 


3.1 Morphological characteristics of A. cerana 
cerana in the Yimeng Mountains 

The bees were dissected and were made 
specimens. We measured and counted the data. The 
results showed that the body length of A. cerana 
cerana of seven populations from the Yimeng 
Mountains was 12.064 —13.351 mm, which was higher 
than the reported value in the literature (11. 0 — 
13.0 mm) (China National Commission of Animal 
Genetic Resources, 2011). The average body length 
of the XL-DLZ population was 13. 351 mm. The 
average forewing length of the seven populations was 
8.198 — 8. 694 mm and the average forewing width 
was 2. 828 — 2. 984 mm. The average cubital index 
was 3. 239 – 3. 688. АП the morphological information 
was presented in Table 4. 
3.2 Morphological index discriminant analysis 
of A. cerana cerana populations in the Yimeng 
Mountains 

The Wilk’ s lambda value from the stepwise 
discriminant analysis was used to classify the 36 
morphological indexes from the 114 broods. Then, a 
discrimination function scatter diagram was drawn 
(Fig. 1). A. cerana cerana populations from different 
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positions were clustered based on the centroid of the 
different populations to estimate the morphological 
distances. As shown in the Fig. 1, the BHY and 
DTZ-ZJZ populations clustered together, while the 
XL-DJW,  XL-DLZ and  SI-AQ populations 
clustered. However, the XZ-XLZ population was far 
from the other populations, and the different broods 
in this population were spread out. 
3.3  Intra-population genetic variation of A. 
cerana cerana populations in the Yimeng 
Mountains 

A total of 58 alleles were detected in the seven 
populations in the Yimeng Mountains with 11 
microsatellite markers. The number of alleles, size 
of the alleles, observed heterozygosity ( Ho) and 
expected heterozygosity ( He) and PIC for each locus 
were listed in Table 5. 

Table 5 Genetic variation information in seven Apis 


cerana cerana populations of the Yimeng Mountains 
based on 11 microsatellite markers 





Epcus Qe x dide (bs ) He Ho PIE 
P 14 158-172 0.8456 0.9115 0.8233 
y 17 127-152 0.7533 0.7788 0.7202 
X 3 111-113 0.2717 0.0000 0.2384 
5 2 147-148 0.0519 0.0000 0.0504 
T 4 132-135 0.5784 0.5841 0.4853 
wW 2 264-265 0.0175 0.0000 0.0172 
K 2 257-258 0.0346 0.0000 0.0339 
L 2 151-152 0.0175 0.0000 0.0172 
M 6 140-156 0.4337 0.3009 0.4072 
R 3 142-148 0.3949 0.4771 0.3395 
Q 3 102-106 0.0270 0.0273 0.0265 
Mea 5.2727 0.3115 0.2800 0.2872 


He; Expected heterozygosity; Ho: Observe heterozygosity; РІС: 


Polymorphism information content. For locus infromation, see Table 2. 


The number of alleles per locus ranged from 2 
to 17, and the average number of alleles observed in 
the 11 microsatellite loci was 5.2727. Across the 
populations, both the W and L loci had the lowest 
Ho, He and PIC. In contrast, the P locus had the 
highest Ho, He, and PIC values (0. 8456, 0.9115 
and 0. 8233, respectively). The average value of 
Ho, He and PIC were 0.3115, 0.2800 and 0. 2872, 
respectively. 

The average number of alleles and Ho and He 
for each population across the 11 loci were shown in 
Table 6. The BHY population had the highest 
average number of alleles. Across the 11 loci, the 
XL-DLZ population had the highest Ho value 
(0.3119), and BHY had the highest He value 
(0.3397). It was interesting that the ST-AQ 


population had the lowest average number of alleles, 
Ho and He. 


Table 6 Average genetic variation information in seven 
Apis cerana cerana populations in 
the Yimeng Mountains 





Populations Number of alleles Average Ho Average He 
XL-DLZ 2.8182 0.3119 0.2856 
XL-DJW 2.9091 0.2611 0.2715 

ST-AQ 2.4545 0.1798 0.1916 
ZZ-ZWZ 3.0000 0.2536 0.2790 
DTZ-ZJZ 2.7273 0.2603 0.2688 
XZ-XLZ 2.9091 0. 2673 0.2781 

BHY 4.0000 0.2823 0.3397 


3.4  Inter-population genetic diversity of A. 
cerana cerana populations in the Yimeng 
Mountains 

The results from the F-statistics analysis were 
shown in Table 7. The genetic diversity of seven 
populations was tested by using the Fit, Fst and Fis 
values of each locus. The negative Fis values for 
some populations indicated an excess of heterozygous 
genotypes with respect to the expected value. For the 
whole population, the inbreeding index ( Fst) was 
0.048 ( P «0. 05); all of the loci supported this 
result, except for the W and L loci. The average Fit 
value was 0. 082. 


Table 7  F-statistics analysis of seven Apis cerana 
cerana populations in the Yimeng Mountains 





Locus Fit =F Fst = 0 Fis =f 
P -0.079 *** 0.041* -0.124 *** 
y -0.013 0.056 * -0.074 *** 
X 1.000 *** 0. 069 * 1.000 *** 
5 1.000 *** -0.027* 1.000 *** 
T -0.045 *** 0.045 * — 0. 094 *** 
W 1.000 *** 0.004 1.000 *** 
K 1.000 *** -0.024* 1.000 *** 
L 1.000 *** 0.004 1.000 *** 
M 0.294 *** 0.032* 0.270 ** 
R -0. 163 *** 0.071 * -0.252“ 
0 -0.006 *** 0.027 * - 0.034 *** 

Mean 0.082 *** 0.048 * 0.0364“ 

Standard deviation 0.114 0.005 0.119 


F The deviation between the actual frequency and theoretical frequency 
of genotype in the total population; 0: The relationship between two 
gametes in each subpopulations; f: The deviation between the actual 
frequency and theoretical frequency in the subpopulation. Mean estimates 
from jack-knife over loci. Significance of F-statistics was done using 
Bonferroni permutations based on 1 000 resamplings, * Р < 0. 05; 
" P«0.01; "* P «0.001. For locus information, see Table 2. 


The Nei' s genetic distance between all the 
populations was calculated according to the allelic 
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numbers and the frequency of the 11 loci (Table 8). 
The Nei’ s genetic distance ranged from 0. 0092 
( XL-DLZ to XZ-XLZ) to 0. 1000 ( XL-DLZ to XL- 


DJW). The genetic identity ranged from 0. 9521 
(BHY to ZZ-ZHZ) to 0. 9908 ( XZ-XLZ to XL- 
DLZ). 


Table 8  Nei' s genetic distance (lower triangle) and the genetic identity (upper triangle) between 
Apis cerana cerana populations in the Yimeng Mountains 





Populations XL-DLZ XL-DJW ST-AQ ZZ-ZHZ DTZ-ZJZ XZ-XLZ BHY 
XL-DLZ 0.9901 0.9805 0.9715 0.9724 0. 9908 0.9790 
XL-DJW 0. 1000 0.9725 0.9757 0.9707 0.9898 0. 9766 

ST-AQ 0.0146 0.0279 0.9636 0.9603 0.9776 0.9532 
ZZ-ZHZ 0.0289 0.0246 0.0371 0. 9588 0.9691 0.9521 
DTZ-ZJZ 0.0280 0.0298 0.0405 0.0421 0.9775 0.9738 
XZ-XLZ 0.0092 0.0102 0.0227 0.0314 0.0227 0.9803 

ВНҮ 0.0212 0.0237 0. 0479 0. 0491 0. 0266 0. 0199 


3.5 Cluster analysis of А. сеғапа cerana 
populations in the Yimeng Mountains 

In this study, a dendrogram was created using 
method to show the 
phylogenetic relationships among the 


the — neighbour-joining 
seven А. 
cerana cerana populations based on Nei' s genetic 
distances (Fig. 3). 

The genetic diversity of the seven populations 
was estimated with the Structure program. The 
length of the burn-in period was 20 000, and the 
number of MCMC reps after burn-in was 100 000. 
The predicted clusters (К value) were 2-6, and 
the number of iterations was 10. The result was 
submitted to Structure Harvester to obtain the 
clusters ( optimum K value). In this study, the 
cluster was K 24 (Fig. 2). The population genetic 
structure figure was obtained with the Disturct 
program (Fig. 3). As shown in the figure, the XL- 
DLZ, ST-AQ and XZ-XLZ populations were 
classified into the same cluster. The DTZ-ZJZ and 
BHY populations were also grouped. Additionally, 
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Fig. 2 The predicted clusters of Apis cerana cerana 
populations in the Yimeng Mountains by 
Structure Harvester 
The figure was obtained from online tool Structure Harvester (http; // 
taylor0. biology. ucla. edu/struct_harvest/). К represents the number 
of population clusters, and the Delta K was the result of software. 


the XL-DJW and ZZ-ZHZ populations were classified 
into the same cluster. 


4 DISCUSSION 


The Yi and Meng Mountains were used as 
coordinates for the Yimeng Mountains. In this study, 
the seven sampling localities are mainly located 
northeast of the Meng Mountain ridge and the Yiyuan 
The altitude of these 
localities ranged from 170 m to 550 m, and these 
localities covered almost all the geographical features 
in the Yimeng Mountains. 

4. 1  Morphological specificity of A. 
cerana in the Yimeng Mountains 

In this study, A. cerana cerana in the Yimeng 
Mountains had displayed longer body lengths, 
shorter tibias, a broader basitarsus and a shorter 


country mountain basin. 


cerana 


forewing. All these morphological characteristics 
were different from documented records ( China 
National Commission of Animal Genetic Resources, 
2011). We deduced that the reasons for the 
morphological related to the 
geographical features and ecological environment in 


specificity were 
the Yimeng Mountains. The Yimeng Mountains 
belongs to a temperate continental monsoon climate 
with abundant nectar and pollen plant resources. 
Studies have shown that a broader basitarsus was 
beneficial for collecting pollen ( Hao, 2014 ). 
Additionally, with the altitude and 
increasing, the body colour of A. cerana cerana was 
lighter and lighter ( Tan et al., 2003; Wang, 
2007). The longer body length was related to the 
high altitude, latitude апа 
temperature. These results correspond with the 
results of Ruttner et al. (2000). 

4. 2 Morphological diversity of A. 
cerana in the Yimeng Mountains 


rainfall 


annual mean 


cerana 


The discriminant analysis results showed that 
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DTZ-ZJZ 
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XL-DJW 
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Fig. 3 Dendogram of the phylogenetic relationships and population structure of seven Apis cerana cerana 


populations in the Yimeng Mountains 


The phylogentic relationships is based on Nei’ s genetic distance using neighbor-joining method. The test of phylogeny was carried out using a bootstrap 


analysis of 1 000 replications; the numbers at the nodes are percentage bootstrap values. 


the seven populations could be classified into three 
classifications. The BHY and DTZ-ZJZ populations 
clustered together. We thought this was because the 
two populations are located northeast of the Meng 
Mountain and the distance between the localities was 
short. Тһе XZ-XLZ_ population 
differences from other populations, and different 


showed great 


colonies of its population also had significant 
differences. The XL-DJW, XL-DLZ, ST-AQ and 
ZZ-ZHZ populations were all located in the basin of 
Perhaps the 
morphological characteristics were derived from the 


the Yimeng Mountains. similar 
similar geographical conditions. 

4.3 Genetic diversity and relationship analysis 
with microsatellite markers 

We selected 11 microsatellite markers from 
more than 3 000 microsatellite markers based on Ji et 
al. (2014) transcriptome sequencing results. The 
11 microsatellite markers used in the present study 
are randomly distributed across different chromosomes 
or linkage groups in the honeybee genome to 
maintain comparability and representativeness. The 
alleles with the highest frequency were the most 
conservative alleles in the populations, and the other 
alleles were produced by gene mutations during 
evolution. Therefore, microsatellite polymorphisms 
could indicate the evolutionary history of the 
species. 

The PIC value is a good measure of gene 
fragment polymorphisms, with the following ranges: 
(1) PIC » 0.5, the locus is a highly polymorphic 
locus; (2) 0. 25 < PIC < 0. 5, the locus is a 
medium polymorphic locus; and (3) PIC <0. 25, 


the locus is a low polymorphic locus ( Vanhala et al., 
2001). The PIC value is related to the availability 
and utilization efficiency of a marker; the higher the 
PIC value of a marker, the higher the heterozygote 
frequency in one population, and the more genetic 
information it provides. In this study, two loci ( P 
and у loci ) 
exhibited high polymorphic characteristics, while 
three loci (T, M and R loci) showed a medium level 


of polymorphism and six loci (X, S, W, K, L and 


among the 11 microsatellite loci 


Q loci) showed a low level of polymorphism ( Table 
5). The mean PIC value across all the loci was 
0.2872 ( Table 5), 


information for the assessment of genetic diversity. 


which can provide some 


The effective number of alleles is also a good 
measure of the genetic variation, especially in 
conservation genetics studies. The average number of 
the alleles was 5. 2727 across the 11 microsatellite 
loci in the present study. However, the number of 
alleles in some populations was small. Still, this 
indicated that the polymorphism 
provided by these 11 
microsatellite loci in the seven populations was rich. 


result also 
information content 

Gene heterozygosity, also called gene diversity , 
is a suitable parameter for investigating genetic 
(1993 ) defined that a 


polymorphic locus must have a heterozygosity of at 


variation. Neuman 
least 0.10. However, in our study, five loci among 
the 11 microsatellite loci exhibited low gene 
heterozygosity. The result suggests that we may have 
selected 11 microsatellite markers that were not 


cerana cerana in the Yimeng 


highly 


appropriate for A. 


Mountains, ав they меге conserved 
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hereditarily. During the long period of evolution, the 
bees from the Yimeng Mountains hardly mutated 
because of the geographic isolation. The mean He 
value can approximately reflect the genetic structure 
variation; for example, the BHY population had the 
highest genetic variability (0. 3397). 

In our study, the average genetic differentiation 
(Fst) value among the broods was 0. 048, which 
indicated that there is no great differentiation among 
the populations. This result strongly supported the 
close relationships among the populations. The 
coefficient Fis value indicates the degree of departure 
from random mating. Positive Fis values indicate a 
significant heterozygote deficit, while negative Fis 
values indicate an excess of heterozygote genotypes 
with respect to the expected value. In our study, the 
average Fis value was 0.036. However, the values 
of five loci among the 11 microsatellite markers were 
1.000, which also indicated the close relationship 
among the populations. Additionally, some loci 
showed a significant deficit of heterozygotes. 


The genetic distance could reflect the 
evolutionary relationship between | different 
populations, and its calculation methods аге 


multifarious ( Jiang et al, 2007; Yang et al., 
2007). The Nei’ s and Reyonld' s genetic distance 
calculation methods are usually used to study 
population genetic variation using microsatellite DNA 
polymorphisms. In this study, the Nei’ s genetic 
distance between all the populations was calculated. 
The Nei’ s genetic distance ranged from 0. 0092 
( XL-DLZ to XZ-XLZ) to 0. 1000 ( XL-DLZ to XL- 
DJW) , and the highest value of genetic identity was 
0. 9908 ( XZ-XLZ to XL-DLZ). The results showed 
that the genetic relationship between XL-DLZ and 
XL-DJW was 


geographic distance between the two populations is 


close. We speculate that the 
small, leading to high gene flow. 

We constructed a phylogenetic tree based on the 
Nei' s genetic distances. The genetic diversity of the 
seven populations was also estimated with the 
Structure program. The results showed that all the 
populations could be grouped into three categories. 
Interestingly, even though the geographical distance 
between the XL-DLZ and XL-DJW populations was 


close, the genetic distance between the two 
populations was far. We believe that the source of 
the two populations is different despite of the high 
gene flow. 

The results of genetic diversity were slightly 
different from the results of morphological index 


Although the XL-DLZ and 5Т-АО 
different from the ХІ-ХІХ 


analysis. 
populations were 


population in morphological analysis, the XL-DLZ, 
the ST-AQ and the XZ-XLZ populations were in the 
same category in microsatellite analysis. In fact, it is 
phenotype was 
determined by genotype. In our study, the clustering 


generally ^ accepted that the 
results according to genetic distance were not in 
consistence with that of the analysis of morphology. 
We thought the reasons were that the samples were 
not enough or the choice of microsatellite loci was 
incomplete. So, the RNA-seq and the choice of 
microsatellite loci would be our next work. 
4.4 Conclusion 

On the whole, the 


morphological and microsatellite marker analysis 


results from the 


showed that the genetic relationship between A. 


cerana cerana populations from the Yimeng 
Compared to the previous 
study 


characteristics 


Mountains was close. 


research results, our revealed different 


morphological and microsatellite 
analysis for A. cerana cerana from the Yimeng 
Mountains. The special geographical environment of 
the Yimeng Mountains and the long history of 
breeding contributed to the diversity and specificity. 
Thus, these A. 


unique strain characteristics and should have higher 


cerana cerana populations have 


levels of protection and utilization as compared to 
other populations. 
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